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ABSTRACT 

The North Polar Spur (NPS) is the brightest filament of Loop I, a large 
circular feature in the radio continuum sky. In this paper, a model consisting of 
two synchrotron emitting shells is presented that reproduces large-scale structures 
revealed by recent polarization surveys. The polarized emission of the NPS is 
reproduced by one of these shells. The other shell, which passes close to the 
Sun, gives rise to polarized emission towards the Galactic poles. It is proposed 
that X-ray emission seen towards the NPS is produced by interaction of the two 
shells. Two OB-associations coincide with the centers of the shells. A formation 
scenario of the Loop I region is suggested. 



Subject headings: ISM: magnetic fields 
neighborhood 



ISM: structure — polarization — solar 



Introduction 



Total intensity surveys reveal a number of radio spurs that can be joined into small 
circles on the sky, so-called radio loops. One of these is Loop I, which has an intriguing 
filament called the North- Polar Spur (NPS). It has been concluded by several authors (e.g. 
Berkhuijsen et al.lll97ll ; lHeileslll979l ; ISalterlll983l . and references therein) that the radio loops 
are correlated with expanding gas and dust shells, energized by supernovae or stellar winds. 
Accordingly, the Loop I superbubble was attributed to stellar winds from the SCO-CEN OB 
association and supernova activity in the same vicinity, with the NPS being the brightest 
segment o f a su pern ova remnant (SNR). Its magnetic field (B-field) has been modelled by 
Spoelstral (119721 ) and iHeiled (119981 ) based on radio and optical polarization data, su ggesting 
that the local B-field is deformed by an expanding shell. According to a model by IWeaver 
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( 119791 ). an SNR, produced by a member of the SCO-CEN association, has expanded inside 
the Loop I bubble. Its shell is just beginning to encounter and interact with the surface 
of the surrounding H I shell (GS 331+14-15). There has been a debate over the expansion 

2 kms -1 , based on data from the 



velocity (v ex p) of GS 331+ 14-15: Weav er gi ves t> exp ~ 



northern hemisphere, whil e iHe iles (19 84) an d lSofue et al.l (119741 ) find a higher t> exp of 19 and 



25 kms 1 . In a picture by iBochkarevi (119871 ) the H I shell has already passed the Sun. 



The relatively low expansion velocity of the neutral gas surrounding Loop I implies an 
age of some 10 6 years and hence suggests the NPS to be part of an old SNR. At this age, 
however, an SNR is well beyond its radiative phase and its optical and radio emission become 
invisible. On the other hand, X-ray and radio emission have been observed t owards the NPS, 
i ndica ting an age at least one order of magnitude less. Based on X-ray data, iBorken fc Iwan 
( 1977 ) suggest that Loop I is an old but recently reheated SNR. This was refined by lEggei 
(119951 ) who argue that Loop I is caused by a shock from a recent supernova, heating the 
inner wall of the SCO-CEN supershell, and giving rise to the prominent X-ray feature of the 
NPS. 



Recently, I Wolleben et al.l (120061 ) published the new DRAO Low- Resolution Polarization 
Survey of the north ern sky at 1.4 GIL 3, which is part of the Canadian Galactic Plane Survey 
(ITaylor et al.ll2003l ). Making use of these new data as well as the WMAP polarization data 
at 23 GHz (Page et al. 2006), a new and detailed study of the Loop I region can be 
performed. In this paper I describe a model consisting of two interweaving magnetic shells, 
which explains some newly detected polarization and depolarization structures, correctly 
predicts the shape of the NPS, as well as providing a clue to the origin of X-ray emission 
observed towards the NPS. 



2. Large Polarized Structures In The DRAO Survey 

The polarization data used in this paper are shown in Figures [T] and [2] The NPS has 
long been known to be one of the most intensely polarized features of the sky. Its filaments 
cover large parts of the northern Galactic hemisphere at latitudes of b > 30°, with an abrupt 
drop of polarized emission below 30° latitude at 1.4 GHz. Surprisingly, structures of the 
NPS are very evident in polarized intensity but its filamentary structure is not detectable 
in polarization angle. In the following, previously unknown polarization features detected in 



1 See also IWollebenl (|2005l ). The data are publicly available and can be downloaded at: 
http : //www. drao . nrc . ca/26msurvey, http://www.mpifr-bonn.mpg.de/div/konti/26msurvey, or from 
CDS, Strasbourg (| Wolleben et al.ll2005h 
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the DRAO polarization survey, which are of relevance for this study, are described: 



2.1. High Latitude Polarized Emission (HLPE) 

In Fig. CD (top), towards the Galactic poles (b = ±90°), a systematic pattern of polarized 
e mission and po l arizat ion angle is visible (HLPE, hereafter). At 408 MHz, the all-sky map 



of lHaslam et al.l (119821 ) shows a weak counterpart of the HLPE in total intensity (Stokes-I) 



with excess emission at the Galactic poles of ~ 4 K (northern p ole) and ~ 5 K (southern 



pole). The Stockert survey at 1420 MHz (IReich fc Reichl 1 19861 ) also shows weak Stokes-I 



emission at the poles of ~ 60 mK (northern pole) and ~ 200 mK (southern pole). These 
temperatures imply a spectral index of —3.4 and —2.6 for the northern and southern HLPE, 
respectively, and thereby suggest that the HLPE is synchrotron emission. 

What could cause enhanced synchrotron emission around the Galactic poles - the 
HLPE? 

First, the HLPE could be (polarized) synchrotron emission from the local spiral arm, 
which would produce a similar pattern of polarization angles towards the Galactic poles. 
Second, the HLPE could be caused by the Galactic halo. However, neither can explain 
the distinct boundaries of the HLPE, which are observed in the second and third Galactic 
quadrant (cf. thin white lines in Fig. [I] top). In fact, synchrotron emission from the halo 
or local spiral arm should be observable at all latitudes where the local Galactic B-field is 
perpendicular to the line-of-sight, and not just at the poles. 



2.2. New Radio Loop 

A previously unknown filament of polarized emission in the southern Galactic hemi- 
sphere (hereafter referred to as "New Loop") is revealed by the DRAO polarization survey. 
The New Loop has an excess polarized intensity of ~ 170 mK at 1.4 GHz relative to its 
surroundings. The elongated shape of the New Loop can be fitted by a small circle. From a 
fit by eye, the center of this circle is at / = 345°, 6 = 0°, with a radius of 65° (dashed line 
in Fig. [T]top). Faint counterparts can be found in total intensity at 408 MHz (~ 5 K) and 
1420 MHz (~ 250 mK), which gives a spectral index of —2.4. This filament can also be seen 
in the WMAP polarization data at 23 GHz. 
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2.3. Depolarization Band 

A minimum of polarized emission from the first Galactic quadrant from ~ 0° to ~ 70° at 
1.4 GHz, within a band confined by remarkably sharp boundaries at b ~ ±30° (cf. Fig[T]top), 
is seen with great clarity for the first time. The percentage polarization within this area is 
around 3%, compared to 20 - 30% outside the band. Obviously, strong depolarization within 
the Galactic plane must "destroy" the polarized emission from the first Galactic quadrant. 
The Depolarization Band is prominent in the DRAO data at 1.4 GHz but not seen in the 
WMAP data at 23 GHz, which ephasises the assumption that it is caused by depolarization 
due to Faraday rotation. The length of the line-of-sight makes differential Faraday rotation 
a likely depolarization mechanisno However, the absence of polarized emission from the 
NPS below b = 30°, in the DRAO data at 1.4 GHz as well as in the WMAP data at 23 GHz, 
cannot be attributed solely to the Depolarization band. It must reflect an inherent property 
of the emission region. This is explained by the model presented in this paper. 



The Model 



The model consists of two synchrotron emitting shells: SI and S2. The shells are 
spherical with constant shell thickness. No emission is produced outside the shells. The Sun 
resides within SI between its inner an d outer surface (see Fig. [3] and JE) ■ T heoretical models 



of magnetic fields in supershells (e.g. iFerriere et al.lll99ll ; iTomisakal Il992l ) predict that the 



ambient B-field is pushed by the shock wave and compressed within the expanding shell. 
The model described here resembles these theoretical B-fields in the simplest way: Magnetic 
field lines run from the polar caps (the magnetic poles) along longitudes of the shell. The 
appearance of the B-field of the shell, that is, its B± and B» components, depends on the 
vantage point. In case of a large nearby shell the projection can result in a complicated 
B-field pattern on the sky. 

Each shell is described by 8 parameters (Tab. 1): the center coordinates I, b, d; the 
inner (r; n ) and outer (r out ) radius of the shell; the angle between the B-field and the line- 
of-sight to the Galactic center (Bg) and to Galactic north {B$)\ and two scaling factors 
describing the intrinsic brightness of each shell. For simplicity, it is assumed that there is 
no Faraday rotation within the shells. The model described here is an attempt to reproduce 
the polarized emission - the emission tracing the regular B-field - because this model does 



2 Differential Faraday rotation occurs if synchrotron emitting and Faraday rotating regions are mixed. In 
this case, polarization vectors with different orientation may superimpose and cancel each other (see also 
Burn 1966). 
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not take irregular or turbulent B-field components into account. Therefore, the model must 
be fitted to polarization maps, as is done here. Moreover, polarization maps at frequencies 
around 1 GHz or less are believed to show rather local emission 

due to depolarization, which helps avoiding confusion with unrelated background emis- 
sion. 

The model uses different field orientations in the two shells. Where the shells overlap in 
space this is obviously not correct. But the complex shape of the edge of the Local Bubble (cf. 
Fig [3]) suggests a very complex evolution, probably the outcome of many earlier stellar winds 
and supernova events. Any previous large-scale magnetic field has probably been tangled by 
this complex evolution. Presumably the Local Bubble is not a unique circumstance, so the 
implication is that most regions in the Galaxy are not characterized by a large-scale quasi- 
uniform component. The two-field model proposed is a simple approximation to a complex 
reality, and is unlikely to be correct over the full extent of the shells. However, the good fit 
to the data demonstrates that it is an adequate model for most of the Loop I region. 

The model was used to calculate the integrated Stokes U and Q values for each line- 
of-sight through the shell complex. Polarized intensity and polarization angle were derived 
from these integrated values, thereby accounting for depolarization due to superposition of 
differently oriented polarization vectors along the line-of-sight (see Appendix for a more 
detailed description). In the northern hemisphere, DRAO data were preferred over WMAP 
data because of the better sensitivity. Some regions in the DRAO maps were masked out 
however: the "Fan-Region" (see Fig. Q}, the two H Il-regions S27 (I = 4°, b = 22°) and S7 
(Z = 350°, b = 22°), and a ±3° strip along the Galactic plane. In the southern hemisphere, 
23 GHz WMAP data (scaled to 1.4 GHz using a spectral index of —3) were preferred over 
the 1.4 GHz polarization survey of Testori et al. (2004) because the western caps lie at b ~ 0° 
where depolarization at 1.4 GHz is likely to be strong. In the WMAP data the Galactic plane 
was masked out because a simple "scaling up" of WMAP data from 23 GHz to 1.4 GHz, 
without taking depolarization effects into account, would result in wrong polarized intensities 
in regions where Faraday rotation along the line-of-sight is high. The southern Galactic pole 
in the WMAP was masked out because here sensitivity is too low to permit accurate fitting. 

Fitting was done by computer, applying an algorithm that searched for the minimum 
of the square-root of the sum of the differences between modelled and observed polarized 
intensities. The algorithm randomly modified model parameters until a good fit was achieved. 
Different initial start values for the fit were chosen to evaluate the uniqueness of the best-fit. 
In order to determine the confidence ranges of the best-fit parameters, models with slightly 
different sets of parameters were tried and visually inspected. 
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4. Discussion 

Figures [2] and [5] show polarized intensity and polarization angle maps of the best-fit 
model. At high Galactic latitudes the HLPE is correctly reproduced by SI. Only SI produces 
HLPE because the Sun is located inside this shell, leading to local emission around the 
Sun. S2 produces the polarized emission of the NPS. At intermediate latitudes, where the 
line-of-sight through the B± component of SI is longest, the New Loop is seen. At low 
latitudes (|6| < 30°) synchrotron emission from the two shells is reduced because the path 
lengths through the shells are short and B± is small. Although the model was only fitted to 
polarized intensities, the polarization angles at intermediate and high Galactic latitudes are 
remarkably well reproduced. 

The predicted western cap region does not fit the observations as well as the eastern cap 
region. However, the western caps are likely to be inside the Local Bubble, while the eastern 
caps may have expanded into a denser medium. Hence, the "real" shells are not likely to 
be of perfectly spherical shape. However, a comprehensive model which includes the full 
complexity of the ISM around Loop I is beyond the scope if this paper, and is probably 
impossible for want of adequate data. 

Portions of the two proposed synchrotron shells currently overlap in space, which, pro- 
jected onto the sky, results in a ring-like region in the northern Galactic hemisphere (see 
Fig. [5] bottom). This region roughly agrees with the location of 1.5 keV X-ray emission from 
the NPS. It is therefore suggested that the X-ray emission is produced by recent interaction 
of the two shells within this region. 



5. Possible Formation History 

The SCO-CEN association can be divided into three subgroups: Lower Centauru s Crux 



(LCC), Upper Centaurus Lupus (UCL), and Upper Scorpius (US). iMaiz-ApellanizI (120011 ) 
calculated the positions of the center of each subgroup in the past, taking the effects of solar 
motion, Galactic rotation, and motions in the z-direction into account. Accordingly, the 
stars whose paths come closest to the center of SI are those of the LCcfi which crossed this 
point 6 ± 2 Myr ago at a distance of 70 pc from the Sun. Closest to the center of S2 is the 
US group at its current location (I = 350°, b = 20°, 145 pc away from the Sun). Stellar 
activity within the LCC subgroup started 11 to 12 Myr ago, and 5 to 6 Myr ago within the 



3 Note that the scattering of members of these subgroups across the sky means that the UCL is almost 
as good a candidate. 
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US fide Geus et al.lll98a ) 



The picture that emerges from these data can be summarized as follows. The first 
supernovae are expected to take place about 3-5 Myr after formation of an association. 
Thus, about 7 to 8 Myr ago, supernova explosions began to occur within the LCC and, 
subsequently, delivered enough energy to inflate a bubble around LCC - the Sl-shell. The 
US, instead, has only recently reached its evolutionary time scale and has just begun inflating 
the Loop I bubble (S2). The shock front of Loop I (S2) hit the LCC bubble (SI) just recently 
(10 4 years ago or less), giving rise to the X-ray emission observed toward the NPS. 

In this picture SI is about 6 Myr old, based on the positional coincidence with the LCC 
at this time. S2 is 1-2 Myr old, based on the the age of the US. This suggests that SI is more 
evolved than S2, which may explain why almost no observational tracers of SI can be found 
in total intensity surveys. SI seems to be an almost dissolved shell, whose last observable 
remnants are the HLPE and the New Radio Loop. 

The scen ario proposed here has compell i ng; sim ilarities to the models of H I shells in 
this region by IWeaverl (119791 ) and iBochkarev! (119871 ). mentioned earlier in the Introduction. 
The ages derived above and best-fit radii of the model give expansion velocities of 8 and 
32 kms" 1 for SI and S2, respectively, usin g the standard mod el for the kinematic age of 
stellar wind bubbles with t^ in = 0.6R/v exp ( 1 Weaver et al.lll977l ). Taking into acc ount the 
uncer tainties in this esti mate, these yelocities agree w ith the velocities of 2 kms -1 ( jWeaver 
19791 ) and 19-25 kms -1 (lHeileslll984l ; ISofue et al.lll974l ) found for H I gas in this region. This 
may imply that the H I distribution in this region is the result of two superimposed H I 
shells that are expanding with different velocities. 



6. Conclusions 

Based on an analysis of the new DRAO polarization survey and scaled WMAP data, a 
model consisting of two synchrotron emitting shells is developed that reproduces large-scale 
structures in the polarized sky. One of these shells, SI, has reached the Sun and gives rise 
to polarized synchrotron emission at high Galactic latitudes, the HLPE described in this 
paper. Where the path length through SI is longest and its B-field is perpendicular to the 
line-of-sight over its whole length in the emitting region, a new radio loop is detected in 
polarized intensity. A younger shell, S2, correctly reproduces emission from the NPS. A 
scenario is proposed in which S2 has recently started to interact with SI, possibly giving 
rise to the observed X-ray emission. The model also predicts correctly the low polarization 
from SI and S2 within the Depolarization Band. The picture of the local ISM suggested in 
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this paper is in agreement with previous studies outlined in the Introduction, although the 
two-shell geometry of the model makes a comparison with previous models difficult. 
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A. Magnetic Field Used In The Model 

In this model an approximation for the magnetic field Bi inside an expanding shell is 
used. Its strength and direction depend on the direction (9, ip) of the ambient magnetic field 
B as well as the inner and outer radii of the shell (r in ,r out ). Looking at Fig. El it is easy 
to see that the magnitude of Bi is maximal / zero where the expansion is perpendicular / 
parallel to B. For \x\ > r out and \x\ < r in the magnetic field strength Bi is assumed to be 
zero. 

The direction of Bi is given by 

Bi = A> ( A1 ) 

\Pi\ 

where pi = b — (b ■ x)x with b the unit vector of the ambient field B and x the unit vector 
pointing towards the i-th line-of-sight element. The field strength of Bi can be evaluated to 

rn = (b- x)\ (A2) 

giving rise to the Stokes parameters 

U{ oc rii B^ 1 ^ 2 sin 2 arctan ^ Jj- 



d(7+1)/2 

q.i oc rii By cos 



2 arctan ' 1,1 



(A3) 



for the i-th line-of-sight. Here B± = a/ Bf + B^ and Bi,Bt are the projections of Bi onto 
the line-of-sight. An energy spectral index of the synchrotron radiation of 7 ~ 2.8 is taken. 
Assuming no intrinsic Faraday rotation, the polarized intensity and the polarized position 
angle is finally found 

PI = y/XP + Q 2 , 



PA = I arctan f ^ 



(A4) 
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where U = J2 u i an d Q — ^2 Qi- 
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Table 1. Best-fit parameters of the model. 
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Fig. 1. — Map of polarized intensity in units of mK (top) and polarization angle (bottom) in 
Galactic coordinates, taken from the DRAO polarization survey at 1.4 GHz (northern hemi- 
sphere) and the WMAP polarization data smoothed to 4° resolution (southern hemisphere). 
The data are shown in rectangular projection to expose structures around the Galactic poles. 
The polarized intensity map exhibits three extended features: the NPS (from I = 310° to 
50° and b = 30° to 80°), the New Radio Loop (centered at I = 40°, b = -55°, about 40° 
in diameter), and the so-called "Fan- Region" (centered at I = 140°, 6 = 5°, about 60° in 
diameter). The solid lines mark the observable extent of the HLPE in the second and third 
Galactic quadrant. The dashed line shows a circle fitted through the New Loop and one of 
the two extended "arms" of the NPS (at about b = 70°). Black contours in the bottom panel 
show 1.5 ke V X-ray emission asso ciated with the NPS at levels of 200 and 350 x 10~ 6 counts/s 
(taken from lSnowden et al.lll995l ). 
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Fig. 2. — Observed (left) and modelled (right) maps of the northern (top) and southern 
(bottom) Galactic poles with polarization vectors overlayed. The contours show polarized 
intensity from the DRAO and WMAP data from 50 to 400 mK, in steps of 50 mk. 
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Fig. 3. — The sketch displays cuts through the model viewed looking down on the Galactic 
plane towards negative latitudes (top) and the vertical plane through the Sun and perpen- 
dicular to the line-of-sight looking towards I = 90° (bottom). The two shells are indicated by 
solid (SI) and dashed (S2) lines. The Sun is indicated by the circled dot. Thin lines indicate 
the B-field orientation of each shell (the B-field component parallel to the image plane). The 



stipp led region shows the Na I distribution around the Local Bubble (taken from ISfeir et al. 



19991 ). Filled triangles and squares show the centers of the LCC and US su bgroup today, 
5 Myr ago, and 10 Myr ago (only for LCC) (taken from iMaiz-Apellania 120011 ) . 
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Fig. 4. — An "engineering drawing" to help the reader depict the 3-dimensional structure 
of the two proposed shells. The solid line indicates SI, and the dashed line indicates S2.In 
contrast to Fig. [3] these drawings show projections of the two shells rather than cuts through 
the Galaxy. The Sun is located at the origin of the coordinate system. The figure shows 
the shells as seen from the Galactic anti-center towards the Galactic center (top left), seen 
sideways (top right), and seen from above (bottom). 
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Fig. 5. — Modelled maps of polarized intensity in units of mK (top) and polarization angle 
(bottom). The dashed line from Fig. [1] is repeated in the top panel. The black contour in the 
bottom panel indicates the region where the SI and S2 shells overlap in space. The model 
was merely fitted to polarized intensity which means that "Fig. top" (model) was fitted 
to "Fig. [[] top" (survey). Nevertheless, the predicted polarization angles ( "Fig. [5] bottom" ) 
resemble the observed pattern ( "Fig. [1] bottom" ) remarkably well. The Fan- Region is not 
part of the model described in this paper. 



Fig. 6. — Geometry of the spherical shell. r m and r out are the inner and outer radii of the 
shell, respectively; B is the background magnetic field and b is its unit vector; Bi is the 
magnetic field inside the shell along the i-th line-of-sight element; and x is the unit vector 
pointing towards the i-th. line-of-sight element. Region I and II indicate regions with strong 
and weak B-fields, respectively. 



